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Yuji Kodama (Ohio State University)* 
Abstract 

This talk gives a survey of our recent studies on soliton solutions of the 
Kadomtsev-Petviashvili equation with an emphasis on the Mach reflection 
problem in shallow water. 



1. Shallow water waves: Basic equations 

We consider a surface wave on water which is assumed to be irrotational and incom- 
pressible (see for examples |2H H] ) . We are mainly interested in a long wave phenomena, 
and assume zero surface tension, that is, we ignore the gravity-capillary waves. Let us 
^ denote the following scales: 

^ Aq ~ horizontal length scale = typical wave length 

ho vertical length scale = asymptotic water depth 
O ao ~ nonlinear scale = typical wave amplitude 

The non-dimensional variables {x, y, z, t, rj, (f)} with the corresponding physical variables 
|Ih^ with the space variables (x, y, z), the amplitude of the wave f], and the velocity potential 

^ are defined as 

d 



X = Aqx, y = Xoy, z = h^z, t = — t, 

Co 

V = aoT], (p = — AoCo0. 

ho 

^ Then the shallow water equation in the non-dimesional form is given by 

O (pzz + (^^<P = 0, for < 2 < 1 + a?7, 

2 '^- = 0' at z = 0, 

>• ^ ^ dX z = 1 + ar], 

j>P Vt + «V0 ■ Vr/ = -0^, 

S3 

where V = {-^, ^) and A = V^, the two-dimensional Laplace operator, for (x, y) G 
The parameters a and /3 are given by 

a = ^ and = f 

ho W 

The weak nonlinearity implies a ^ 1, and the weak dispersion (or long wave assump- 
tion) implies /3 ^ 1. With a small parameter e ^ 1, we assume a ~ /3 = (9(e). 
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First we note, from the first two equations, tliat can be written formally in the 
form, 

y, z, t) = cos (z v^) ^ = V - /5y AtA + P^Ya^^^ + 

where ip{x,y,t) = 0(3^, y, 0, t). Then the equations at the surface z = 1 + ar] give a 
Boussinesq-type system, 

+^ ((A^)^ + ■ V(A^) - 2vA^p,) + ^A^^, = Oie'). (1.2) 

rjt + Aip+aV ■ {rjV^lj) - /J-A^^ 

6 

(Vr^ ■ V(A^) + r^A^^) + ^A^^ = 0(6=^). (1.3) 

We now derive the KP equation with higher order corrections which will be a key 
equation for a physical application. The KP equation is obtained under the assumption 
with a weak dependence in the y-direction (quasi-two-dimensionality), and we introduce 
a small parameter 7 so that the y-coordinate is scaled as 

C:=v^l/, with i = 0{e). (1.4) 

We also consider a far field with the scaled coordinates ^ = x — t and t = et (i.e. a 
unidirectional approximation). Then eliminating rj in the equations (1.2) and ( |1.3 ), we 
obtain the KP equation for the function v := ipS.iCXy^) with higher order corrections 
up to C(e2), 

/3 _i 

19 02 o 53 \ /37 72 

+ + «/5 ( + ) + ~ T ^^^^^ 

+ a7 Q^^^'^CC + 2t;c/^-'^;c - Id-\v^)c^ + ^v^D^S^^ = 0{e'). (1.5) 

where := is a formal integral operator. The first line of this equation of the 
order 0{e) is the KP equation, and the terms in the second and third lines are the 
e^-order corrections to the KP equation. The wave amplitude r] is then given by 

V = v + ^v'- ^v^^ + ^D~\^ + 0{e'). (1.6) 



Remark 1.1 In terms of physical coordinates, the KP equation is given by 



3Co _ _ Cg/in _ 

2ho 6 



I '^0 ~ n 



As a particular solution, we have a solitary wave solution in the coordinate perpendic- 
ular to the wave crest, x = x cos \E'o + y sin \E'o; 



fj = aosech 



3ao 



ihl cos2 ^0 



, CLq 1 2 T I 

X - Co COS Wo ( 1 + — + - tan Wo ) t - Xo 



;i-7) 



where oq > 0, \l/o and xq are arbitrary constants. Recall that the KP equation is derived 
under the assumption of quasi-two dimensionality, that is, 7 = tan^ \l/o = 0{e), and 
the solution (1.7) becomes unphysical for the case with a large angle. That is, the width 
of the solitary wave depends on the propagation direction, \l/o, which should be corrected 
when we apply the KP equation to a physical problem. In this talk, I will explain how 
this can be fixed using a normal form theory concerning the higher order corrections to 
the KP equation, and show that some of the exact solutions of the KP equation can be 
used to describe two-dimensional interaction phenomena in shallow water, called the 
Mach reflection. 

2. Normal form for the KP equation with higher order correc- 
tions 



In order to study the higher order corrections, we define a normal form for (1.5). For 
this purpose, we first put (1.5) into a canonical form of the KP equation with the 
change of variables, 

2 3a 



Then (1.5) becomes 

4:Ut = 6uux + Uxxx + 3D~^uyy 

19 5 53 3 3 

+ -^UXXXXX + 7,UUXXX + -^UxUxX + ^UXYY - -D'^UyVYY 

oU 3 o 2 4 

Q 

+ -uD-^Y + 4:UyD-^uy - -D-\u^)yy + uxD-\yy + C(et). (2.1) 

Note here we have the orders u ~ 0{e), dx ~ C'(e^), ~ 0{e^^), and 9y ~ 0{e). 
Here the new variables (X, Y, T) are related to the physical ones with 

x-coi=hoX, y = hoY, i=-^T. (2.2) 

/Co 

The wave amplitude t] in terms of u is given by 

2 12 1 

arj = -u + -u^ - gUxx + -D~'^uyy + 0{e^) 

Hereafter we use the lower case letters {x,y,t) for {X,Y,T), and the KP variables 
can be converted to the physical variables directly through the relations (2.2). 

2.1. The normal form for the KdV equation 

Before discussing a normal form for (1.5), we give a brief summary of the result in 
[HI [ini [S] for the case of the KdV equation. Taking dyu = in (2.1 ), we have the KdV 
equation with higher order corrections, 

/19 5 53 \ 9, , , 

AUt = 6UU:^ + + I —U^xxxx + ^UU^xx + yM^Maiz I + 0[e^). (2.3) 



In [01 [ini E], we found that one can transform (2.3) with a formal change of variable, 
a Lie exponential transformation e^'^ with the generating function ip (see [5]), 

u = e''- oU = U+ (u,, + ^U' + hj,D-^u\ + C(e3), (2.4) 



into the equation which we referred to as the normal form of (2.3) 



19 

4Ut = 6UU, + U,,, + — (f/™ + lOUU^^, + 20U,U,, + 30U^U,) + 0{e 

oU 



The point here is that the higher order term of ^(ea) in this equation is the 5th order 

7 

symmetry of the KdV equation, and hence the normal form is integrable up to ^(es). 
That is, we have an integrabihty not only at the KdV of ^(ea) but also at the next 

7 

order correction of (^(es). This is true even for the general form of the higher order 
correction, 



with arbitrary coefficients ai, . . . , a^. This implies that any weakly nonlinear long-wave 
equation whose leading order is approximated by the KdV equation is asymptotically 
integrable up to the next order approximation [9l [TOl E]. 

Note also that the normal form admits one-soliton solution in the form. 



U = Ao sech^ 




'x + xo{t)), 



where Xo{t) is determined by ^ 
of the higher order KdV equation (|2.3l) is given by the transformation (|2.4l), i.e 



ChKdv = + + 0{e^). Then the solution 



u = AoS' + ( AlS' 



where S := sech 
is given by 



— {x + Xq). Notice that the amplitude t] of (1.6) for the KdV case 



av = -A,S' + -AlS' + -AlS' + 0{e') 



(2.5) 



2.2. Normal form of the KP equation 



Now we define a normal form of (2.1) by a Lie exponential transform, 

u = e''-oU = U+ {PiU^x + /32t/2 + l3,,UxD-^U + P^D-^Uyy) + 0(e=^), 

where /5j's are determined such a way that the transformed equation (normal form) 
has a "good" property. First, we require that the normal form is reduced to the KdV- 



normal form when we have dyU = 0. That is, from (2.4), we have /3i = 1, 



3,/33 



1 

2' 



but /34 remains free. Then we require that the normal form admits a solitary wave in 
the form of one-soliton solution. 



U = An sech^ 




X + xo{y,t)). 



This determines uniquely (34^= ^, and the normal form of (2.1) is given by 



+ — iu 

60 ^^^^^ 



'yy 



+ IWU^xx + 20f/,f/,, + ?>W^Ux) + -U^yy 



- UD-^a 



yy 



lUyD-^Uy 



4 

^-D-\U\y+^-UxD-'Uyy + 0{e 



yyyy 



(2.6) 



Then the phase xo{y,t) satisfies 



Setting the angle \&o of the soliton with ^ = tan^o, we have the velocity ChKP = ^ 
of the KP soliton with the higher order corrections. The corresponding solitary wave 
solution u is then given by 

u = AoS' + (aIS' + ^Aotan^ ^H.S' - IaIS') + 0{e^) 



and the wave amplitude ?7 of (1.6) is given by 

= Ia,S' + \ tan^ ^>,S' + ^AlS' + \aIS' + 0{e^) 

= ^[Ao]S' + ^[AorS' + ^[AofS' + 0{e'), (2.7) 



where [Aq] := Aq{1 + tan^ \E'o) = ^o/cos^\E'o (cf.(1.7)). Notice that the correction to 
quasi-two-dimensional approximation can be absorbed into the amplitude of the KdV 
equation (cf.(2.5)). This formula gives the relation between the observed amplitude rj 
from numerical simulation (or experiment) of shallow water wave system and the KP 
amplitude Aq. 

3. The KP solitons 

Here we give a brief summary of the soliton solutions of the KP equation [21 HH [12], 

- 4:Ut + 6uu^ + u^^^ + D'^^Uyy = 0. (3.1) 

We write the solution of the KP equation in the r-function form, 

u{x,y,t)=2dllnT{x,y,t). (3.2) 

where the r-function is assumed to be the Wronskian determinant with N functions 
/i's (see for examples [221 16]). 

r = Wr(/i,/2,...,/^). (3.3) 

Here the functions {/i, . . . , Jn} form a set of linearly independent solutions of the linear 
equations, 

dyfn = d^fn, dtfn = d^fn- 

In particular, we consider a set of finite dimensional solutions for {/i, /2, . . . , /at}, 

M 

fn{x, y, t) = ^ an,mEm{x, V, t) with Em = 6^™ := exp(K„x + nl^y + K^t). 

m=l 

Thus this type of solution is characterized by the ordered parameters {ni < ■ ■ ■ < km} 
and the N x M matrix A := {an,m) of rank(A) = N, that is, we have 

(/i, /2, . . . , In) = {E,, E2,..., Em)A^ . (3.4) 



Note that {Ei, E2, . . . , Em} gives a basis of and {/i, /2, . . . , /at} spans an A^- 
dimensional subspace of M*^. This means that the A- matrix can be identified as a point 
on the real Grassmann manifold Gr(A^, M) (see [TTj 2J). More precisely, let Mnxm{^) 
be the set of all iV x M matrices of rank N. Then Gr(A^, M) can be expressed as 



Gt{N, M) = GL^(M)\M, 



AfxMl 



where GLAr(]R) is the general linear group of rank A^. This expression means that other 
basis ((?!,..., qn) = (/i, • • • , fN)H for any H G GLAr(]R) spans the same subspace, that 
is, A — 7- A (GLAr(]R) acting from the left). Notice here that the freedom in the 
^-matrix with GLjv(K) can be fixed by expressing A in the reduced row echelon form 
(RREF). 



Now using the Binet-Cauchy Lemma for the determinant, the r-function of (3.3) 
can be expressed in the form, 

T{x,y,t)= J2 ^j{A)Ej{x,y,t), (3.5) 

where Aj(y4) is the N x N minor of the A- matrix with columns marked by J = 
{ji,...,jN} e ('^f'), the set of A^-subindices of the set [M] := {1,...,M}, and 
Ej{x,y,t) is given by 

Ej{x, y,t) = Wt{E,^, ...,E,J = JJlSVn - Sj^ii ■ ■ ■ Ej^. 

l<m 

We are also interested in non-singular solutions. Since the solution is given by m = 
29|(lnr), the non-singular solutions are obtained by imposing the non-negativity con- 
dition on the minors, 

Aj(A)>0, for all Je {^^^ (3.6) 

This condition is sufficient for the non-singularity of the solution for any initial data 
(see [15j for the necessary condition for the regularity). We call a matrix A having the 



condition (3.6) totally non-negative matrix, referred to as TA^A'^ matrix, and the set of 
those matrices forms totally non-negative Grassmannian, denoted by Gr+(A^, M) (see 

e.g. mm)- 

Example 3.1 Let us express one line-soliton solution in our setting. Here we also 
introduce some notations to describe the soliton solutions. One soliton solution is 
obtained by the r-function with M = 2 and N = 1, i.e. r = fi = auEi + ai2-E'2- with 



A = (oii 012). Since the solution u is given by (3.2), one can assume an = 1 and 
denote 012 = a > 0. Then we have 

u = 2dl Inr = ^(ki — ^2)^ sech^ ^(6^1 — ^^2 — In a). 

Thus the solution is localized along the line 61 — 62 = In a, hence we call it line-soliton 
solution. We emphasize here that the line-soliton appears at the boundary of two regions 
where either Ei or E2 is the dominant exponential term, and because of this we also call 
this soliton a [l,2]-soliton solution. We refer to each of these asymptotic line-solitons 



as the [i, j]-soliton. The [i, j]-soliton solution with i < j has the same (local) structure 
as the one-soliton solution, and can be described as follows 



u 



with some constant ©jij]- The amplitude ^[jj], the wave-vector 'K^i^j] and the frequency 
are then expressed by 



Note that K[jj] and ^[ij] satisfy the soliton- dispersion relation, AQ[ij]K^- = {K^-j-^)^ + 
3{K^-j^Y. The direction of the wave-vector 'K.^ij^ = {K^- j-^, K^- j^) is measured in the 
counterclockwise rotation from the y-axis, and it is given by 



tan = + Kj, 



that is, '^[ij] gives the angle between the line K[j ,,•] ■ x = const and the y-axis. 
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Figure 3.1: One line-soliton solution of [z, j]-type and the corresponding chord diagram. 
The upper oriented chord represents the part of [z,j]-sohton for y ^ and the lower 
one for |/ <S 0. 



In Figure 3.1, we illustrate one line-soliton solution of [i,j]-type. In the right panel 
of this figure, we show a chord diagram which represents this soliton solution. Here the 
chord diagram indicates the permutation of the dominant exponential terms Ei and Ej 
in the t -function, that is, with the ordering Ki < Kj, Ei dominates in x 0, while Ej 
dominates in x ^ (see Section^ for the precise definition of the chord diagram). 

4. Classification of soliton solutions 

We now present the classification theorems of the soliton solutions generated by the 



T-function defined in (3.5). Here we consider the matrix A to be in RREF, and we also 



assume that A is irreducible as defined below: 

Definition 4.1 An N x M matrix A is irreducible if each column of A contains at 
least one nonzero element, or each row contains at least one nonzero element other 
than the pivot once A is in RREF. So the irreducibility implies that we consider only 
derangements (i.e. no fixed points) of the permutation. 

Our classification scheme of the soliton solutions is given by identifying the asymp- 
totic line-solitons as y ^ iioo. We denote a line-soliton solution by (A^_, A^+)-soliton 
whose asymptotic form consists of line-solitons as y — )• — oo and line-solitons for 



?/ — )■ oo in the xy-pla.ne. The next Proposition provides a general result characterizing 
the asymptotic line-solitons of the (A^_, A^+)-soliton solutions (the proof can be found 
in [2]): 

Proposition 4.1 Let {ei,e2, ■■■ ,eN} and {gi, g2, ■ ■ ■ , gn-N} denote respectively, the 
pivot and non-pivot indices associated with an irreducible, N x M, TNN A-matrix. 



Then the soliton solution obtained from the T-function in (3.5) with this A-matrix has 
the following structure: 

(a) For y ^ 0, there are N line-solitons of [en,jn]-type for some jn- 

(b) For y '^0, there are (M — A^) line-solitons of [im, gm]-type for some im- 



An important consequence of Proposition AA_ is that it defines the pairing map 
TT : [M] ^ [M] with 

7r(e„) = j„ , n = 1,2, . . . , A^, 

(4.1) 

T^{gm) =«m, m = 1,2, . . . ,M - A^. 

Recall that {en}n=i and {gm}m=i respectively, the pivot and non-pivot indices of 
the y4-matrix and form a disjoint partition of [M]. Then the unique index pairings in 



Proposition 4.1 imply that the map vr is a permutation of M indices, i.e. vr G S'm, the 



symmetric group of permutations. Furthermore, since Ti{en) = j„ > e„, n = 1, . . . , A^ 



and 7r{gm) = im < Qmi m = 1, . . . , M — A^, vr defined by (4.1) is a permutation with 
no fixed point, i.e. derangements. Yet another feature of vr is that it has exactly A^ 
excedances defined as follows: an element / G [M] is an excedance of tt if 7r(/) > /. The 



excedance set of vr in (4.1) is the set of pivot indices {ei,e2,. . . .cn}- Then we have 
the following characterization for the line-soliton solution of the KP equation |2]. 

Theorem 4.2 Let A be an N ^ M , TNN, irreducible matrix which corresponds to a 



point in the non-negative Grassmannian Gr'^(N, M). Then the t -function (3.5) associ- 
ated with this A-matrix generates an (M — A^, N) -soliton solutions. The M asymptotic 
line-solitons associated with each of these solutions can be identified via a pairing map 



TT defined by (4.1). The map vr G Sm is a derangement of the index set [M] with N 



excedances given by the pivot indices {ei, 62, ... , e^} of the A-matrix in RREF. 



Theorem 4.2 provides a unique parametrization of each TNN Grassmannian cell in 



terms of the derangement of Sm- One should note that Theorem 4^ does not give us the 
indices j„ and i^ in the [e„, j„] and [imiOrr] line-solitons. The specific conditions that 
an index pair [i, j] identifies an asymptotic line-soliton are obtained by identifying the 
dominant exponential in each domain in the xy-plane. To visualize those asymptotic 
solitons, we define a chord diagram: 

Definition 4.3 A chord diagram associated to a derangement vr G Sm is defined as 
follows: Consider a line segment with M marked points by the numbers {1, . . . , M} in 
the increasing order from the left. 

(a) If i < TT{i) (expedience), then draw a chord joining i and n{i) on the upper part 
of the line. 



(b) If j > -n{j) (deficiency), then daw a chord joining j and Tr{j) on the lower part 
of the line. 



Then Proposition and Theorem 4^ imply that each chord diagram identifies the 
types of sohtons appearing in the asymptotic regions \y\ ^ 0. It is also useful to 
consider those marked points as the values of /c-parameters, so that for each chord 
joining i and j, its length gives the amplitude and the sum of the joined points gives 
the inclination of the line-soliton of [i, j]-type. Figure 4.1 illustrates the time evolution 
of an example of (3, 3)-soliton solution. The chord diagram shows all asymptotic line- 
solitons for y — ±oo. 
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Figure 4.1: An example of (3, 3)-soliton solution. The permutation of this solution is 
71 = (451263). The /c-parameters are chosen as {ki, K2, ■ ■ ■ , i^e) — ' ^ 1 n 1 1 3 
(See [H] for a classification of those patterns.) 
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5. Shallow water waves: The Mach reflection 

In this last section, we discuss a real application of the exact soliton solutions of the 
KP equation described in the previous sections to the Mach reflection phenomena in 
shallow water. Let us flrst discuss some exact solutions relevant to this phenomena. 

5.1. Exact soliton solutions from Gr+(2,4) 
Recall that the r-function for Gr+(2,4) is given by 

r(x,i/,t)= ^ A,j(v4)Ei,,(x,|/,t) (5.1) 

l<i<j<4 

where Aj ,, (y4) > is the Pliicker coordinates, the 2x2 minors consisting of z-th and 
j-th columns of the 2x4 A-matrix, and i^jj = Wr(e^%e^^) = (k^ — Ki)e^^~^^^ (note 



Eij > with the order Ki < kj). Proposition 4.1 shows that r-function (5.1) generates 
a soliton solution which consists of at most two line-solitons for both y — )• ±00. We 
consider the following two types: one consists of two line-solitons of [1,2] and [3,4] 
for both |F| ^ 0, which is called 0-type soliton ("O" stands for original, see [TT]): 
the other one consists of [1,3] and [3,4] line-solitons for Y ^ and [1,2] and [2,4] 
line-solitons for Y <^ 0. Theorem 14.21 indicates that this soliton can be referred as to 
(3142)-type, because those line-solitons represent a permutation n = (3142)- Figure 



5.1 illustrates the contour plots of 0-type and (3142)-type solutions in the xy-plane. 



and the corresponding chord diagrams. The upper chords represent the asymptotic 
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Figure 5.1: Contour plots and the chord diagrams. Left: 0-type (or (2143)-type) 
solution. Right: (3142)-type solution. The length of [1, 4]-soliton changes in t. 



solitons [i,j] for y ^ 0, and the lower chords for the asymptotic solitons [i,j] for 
y '^0. The A-matrices for those solutions are respectively given by 

. _ A a 0\ . _fl a -c\ 

AO-^Q 1 b)' A3142)- Q ^ ^ J, (5.2) 

where a,b,c> are constants determining the locations of the solitons (see [2]). Notice 
that the r-function for the (3142)-type contains five exponential terms, and the r- 
function for 0-type with c = in (5.2) contains only four terms. 

Let us fix the amplitudes and the angles of the solitons in the positive x regions for 
both O- and (3142)-types, so that those solutions are symmetric with respect to the 
X-axis (see Figure 5.1): 



An 







A[l,2] = 


^[3,4] 


(Otype) 


A[l,3] = 


A[2,4] 


((3142) type) 


-^[1,2] 


= ^[3,4] 


>0 (Otype) 


-^[1,3] 


= ^[2,4] 


>0 ((3142) type 



(5.3) 



Then we express the K-parameters in terms of Ao and tanipo with ki = —k^ and 
ti2 = —t^3 (due to the symmetry): In the case of 0-type, we have 



Kl = 

2 



^ (tan ^0 + V^) , '^2 = - ^ (tan $o - • (5-4) 



The ordering K2 < K3 then implies tanipo > V^A^- On the other hand, for the (3142)- 
type, we have 

= (tan^o + ^2A^) , ^2 = 1 (tan^o - V^Aq^ . (5.5) 

The ordering K2 < K3 implies tan \E'o < V^Aq. Thus, if all the solitons in the positive 
x-region have the same amplitude Aq for both O- and (3142)-types, then an 0-type 
solution arises when tan\l/o > V^A^, and a (3142)-type when tan\l/o < y/2A^. Then 
the limiting value at K2 = K3 (= 0) defines the critical angle \l/c. 



tan^c := V2A^. 



(5.6) 



Note that at the critical angle, i.e. K2 = K3, the r-function has only three exponential 
terms, and this gives a "Y"-shape resonant solution as Miles noted [19]. 

One should note that for (3142)-type solution, the solitons in the negative x-region 
are smaller than those in the positive region, i.e. ^[3,4] = ^[1,2] = | tan^ < Aq, 
and the angles of those in the negative x-regions do not depend on ^0 and ^[3,4] = 
-^[12] = Two sets of three solitons {[1,3], [1,4], [3,4]} and {[2,4], [1,4], [1,2]} are 
both in the soliton resonant state. These properties of the (3142)-type solution are the 
same as those of Miles's asymptotic solution for the Mach reflection in shallow water 
[IP]. However one should note that the [1, 4]-soliton corresponding to the Mach stem 
becomes a soliton solution only in an asymptotic sense. Then the exact solution of 
(3142)-type can provide an estimate of a propagation distance, at which the amplitude 
is sufficiently developed. 

5.2. The Mach reflection in shallow water 

In [19], J. Miles considered an oblique interaction of two line-solitons using 0-type 
solutions. He observed that resonance occurs at the critical angle \l/c; and when the 
initial oblique angle \l/o is smaller than the 0-type solution becomes singular. He 
also noticed a similarity between this resonant interaction and the Mach reflection 
found in shock wave interaction (see for example [HI El]). This may be illustrated by 



the left figure of Figure |5.2[ where an incidence wave shown by the vertical line is 
propagating to the right, and it hits a rigid wall with the angle \l/o. If the angle \l/o is 
large, the refiected wave behind the incidence wave has the same angle \l/o; i-e. a regular 
refiection occurs. However, if the angle is small, then an intermediate wave called the 



Mach stem appears as illustrated in Figure 5^, The Mach stem, the incident wave and 
the refiected wave interact resonantly, and those three waves form a resonant triplet. 



The right panel in Figure 5^ illustrates the wave propagation which is equivalent to 
that in the left panel, if one ignores the effect of viscosity on the wall (i.e. no boundary 
layer). These refiection patterns are well explained by the exact solutions of O- and 
(3142)-types [HI EH]. In the next section, we give a brief history of studies on the Mach 
refiection. 





Figure 5.2: The Mach refiection. The left figure illustrates an incidence wave propagat- 
ing parallel to the wall with the mirror image. The right figure is an equivalent system 
to the left one. The resulting wave pattern shown here is a (3142)-soliton solution. 
The angle $ becomes zero if the initial angle satisfies \l/o > "^c, i-e. no stem. 



5.3. Previous numerical results 

One of the most interesting tilings of the Mach reflection is that the KP theory predicts 
an extraordinary four- fold amplification of the stem wave at the critical angle [I9] . The 
maximum amplitudes for those symmetric cases of 0-type and (3142)-type occur at 
y = (i.e. at the wall), and they are expressed by a parameter k, called the Miles 
parameter defined by (the original parameter was given by /cm = ^o/"\/2yl^), 

tan^o „x 

(5.7) 



V2Ao 

Then from the KP theory, the amplification factor ao := MmaxMo is given by [IHl [2] 

(1 + kf, for A; < 1, 

"°=^ ' for k>l. 



1 + Vl-k 



-2 ■ 



Several laboratory and numerical experiments tried to confirm the formula (5.8), in 



particular, the four-fold amplification at the critical value k = 1 (see for example [20 | [T7 t 
m 123] ) • In 15 ) Funakoshi made a numerical simulation of the Mach reflection problem 



using the system of equations equivalent to the Boussinesq-type system (1.2) and (1.3) 
up to the first order. He considered the initial wave with a small amplitude = 0.05, 
and concluded that his results agree very well with the resonantly interacting solitary 
wave solution predicted by Miles. However his results on the amplification parameter a 
are slightly shifted to the lower values of the original Miles parameter fc^ = ^o/ a/Soj. 
Tanaka in [23] then re-examined Funakoshi's results for higher amplitude incidence 
waves with Oj = 0.3 using the high-order spectral method. He noted that the effect of 
large amplitude tends to prevent the Mach reflection to occur, and all the parameters 
such as the critical angle are shifted toward the values corresponding to the regular 
reflection (i.e. 0-type). For example, he obtained the maximum amplification factor 
a = 2.897 at kM = 0.695. Because of the quasi-two dimensional approximation, i.e. 
|\['o| ^ 1, Miles in his paper [19] replaced tan \E'o by ^E'o, and then in [HES], the authors 
continued on to use this replacement. Then their computations with rather large values 
of gave significant shifts of the parameter k. In the next section, we re-evaluate 
their results using the new parameters obtained from the normal form. 

5.4. Re-evaluation of the previous results using the corrected parameters 
obtained by the normal form 

We recall that the KP equation is derived under the assumptions of quasi-two dimen- 
sionality, weak dispersion and weak nonlinearity. Then for a physical application of 
the KP theory, one needs to include higher order corrections to those assumptions us- 



ing the normal form (2.6). We then convert the observed amplitude in the numerical 



simulations (or the experiments) to the corresponding KP amplitude via the formula 



(2.7). 



Here we just summarize the results: Let be the amplitude of the incident wave 
used in the numerical simulation (or the experiment), and \E'j = \l'o be the inclination of 
the incidence wave. Also let a be the amplification factor obtained from the numerical 
simulation (or the experiment). Then the Miles parameter in terms of the observed 
amplitude is given by 

tanvt/o a/1 + + 5ajtan\E'o ,r 
— — -. (5.9) 



y/2Ao A/6aiCos\E'o 



ao 




Figure 5.3: Numerical results of the amplification factor a versus the parameter k. The 
circles show Funakoshi's result the squares show Tanaka's result [23]. The black 
dots shows the experimental results by Yeh et al [16] . 



The corrected amplification factor (in the KP coordinate) is given by 

a(l + VI + 5ai) 



"0 = ■\ (1 + Vl + 5aai) cos^ vE^q 

a 



if k<\ 
if k>\ 



(5.10) 



We then re-evaluate their results with the new formulae ( 5.9 ) and (5.10), and the results 



are shown in Figure 5.3 Since Funakoshi's simulations are based on small amplitude 
incidence waves, his results agree quite well with the KP predictions. Tanaka's results 
are also in good agreement with the KP theory except for the cases near the critical 
angle (i.e. = 1), where the amplification parameter ao gets close to 3. This region 
clearly violates the assumption of the weak nonlinearity. However, the original plots 
of Tanaka's are significantly improved with those formulae (5.9) and (5.10). The black 



dots in Figure 5.3 indicate the results of recent laboratory experiments done by Yeh 

The laboratory experiments are performed using 7.3 m long 



and his colleagues 

and 3.6 m wide wave tank with a water depth of 6.0 cm. In the talk, I will show the 
details of the experiments including some movies and pictures of several real shallow 
water solitons. I will also discuss a stability problem of solitons solutions, and present 
some numerical simulations which indicate a convergence of KP solutions to certain 



exact soliton solutions similar to the case of the KdV equation [131 [S] (see Figure 5.4). 
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Figure 5.4: Numerical simulation of with the V-shaped initial wave. The solution 
converges to (3142)-type exact solution. Notice a completion of a chord diagram. 
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